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Abstract. New reductions for the multicomponent modified Korteveg-de Vries (MMKdV) 
equations on the symmetric spaces of Dill-type are derived using the approach based on 
' the reduction group introduced by A.V. Mikhailov. The relevant inverse scattering problem 

is studied and reduced to a Ricmann-Hilbert problem. The minimal sets of scattering 
data Tj, i = 1, 2 which allow one to reconstruct uniquely both the scattering matrix and the 
potential of the Lax operator are defined. The effect of the new reductions on the hierarchy 
of Hamiltonian structures of MMKdV and on T, are studied. We illustrate our results by the 
MMKdV equations related to the algebra g ~ so(8) and derive several new MMKdV-type 
equations using group of reductions isomorphic to Z2, Z3, Z4. 

(-h , Key words: multicomponent modified Korteveg-de Vries (MMKdV) equations, reduction 

group, Riemann-Hilbert problem, Hamiltonian structures 
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1 Introduction 

The modified Korteweg-de Vries equation pQ 



O ' It + + 6eq x q z (x, t) = 0, e = ±1, 

OO ' 

has natural multicomponent generalizations related to the symmetric spaces [2]. They can be 
integrated by the ISM using the fact that they allow the following Lax representation: 



d 



= \i— + Q(x, t) - XJj ^(s, t, A) = 0, (1.1) 

Mil) = (ij t + V (x, t) + XVi(x, t) + X 2 V 2 {x, t) - 4A 3 j\ tp(x, t, A) = if>(x, t, A)C(A), 

V 2 {x,t)=4Q(x,t), V 1 {x,t)=2iJQ x + 2JQ 2 , V (x, t) = -Q xx - 2Q 3 , 

where J and Q(x, t) are 2r x 2r matrices: J is a block diagonal and Q(x, t) is a block-off-diagonal 
matrix. The corresponding MMKdV equations take the form 

^ + + 3 (Q X Q 2 + Q 2 Q X ) = 0. 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
[http: / /www.emis.de/journals / SIGMA /symmetry2007.html] 
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The analysis in [2j [31 [3] reveals a number of important results. These include the correspond- 
ing multicomponent generalizations of KdV equations and the generalized Miura transforma- 
tions interrelating them with the generalized MMKdV equations; two of their most important 
reductions as well as their Hamiltonians. 

Our aim in this paper is to explore new types of reductions of the MMKdV equations. 
To this end we make use of the reduction group introduced by Mikhailov [5J [6] which allows 
one to impose algebraic reductions on the coefficients of Q(x,t) which will be automatically 
compatible with the evolution of MMKdV. Similar problems have been analyzed for the iV-wave 
type equations related to the simple Lie algebras of rank 2 and 3 [7J [8] and the multicomponent 
NLS equations [U [10]. Here we illustrate our analysis by the MMKdV equations related to the 
algebras q ~ so(2r) which are linked to the Dill-type symmetric spaces series. Due to the fact 
that the dispersion law for MNLS is proportional to A 2 while for MMKdV it is proportional 
to A 3 the sets of admissible reductions for these two NLEE equations differ substantially. 

In the next Section [2] we give some preliminaries on the scattering theory for L, the reduction 
group and graded Lie algebras. In Section [3] we construct the fundamental analytic solutions 
of L, formulate the corresponding Riemann-Hilbert problem and introduce the minimal sets of 
scattering data Tj , i = 1,2 which define uniquely both the scattering matrix and the solution 
of the MMKdV Q(x,t). Some of these facts have been discussed in more details in [10], others 
had to be modified and extended so that they adequately take into account the peculiarities of 
the Dili type symmetric spaces. In particular we modified the definition of the fundamental 
analytic solution which lead to changes in the formulation of the Riemann-Hilbert problem. 
In Section [J] we show that the ISM can be interpreted as a generalized Fourier [10] transform 
which maps the potential Q{x,t) onto the minimal sets of scattering data Tj. Here we briefly 
outline the hierarchy of Hamiltonian structures for the generic MMKdV equations. The next 
Section [5] contains two classes of nontrivial reductions of the MMKdV equations related to 
the algebra so(8). The first class is performed with automorphisms of so(8) that preserve J; 
the second class uses automorphisms that map J into — J. While the reductions of first type 
can be applied both to MNLS and MMKdV equations, the reductions of second type can be 
applied only to MMKdV equations. Under them 'half of the members of the Hamiltonian 
hierarchy become degenerated [111 [2]. For both classes of reductions we find examples with 
groups of reductions isomorphic to Z2, Z3 and Z4. We also provide the corresponding reduced 
Hamiltonians and symplectic forms and Poisson brackets. At the end of Section [5] we derive the 
effects of these reductions on the scattering matrix and on the minimal sets of scattering data. In 
Section [6] following [3] we analyze the classical r-matrix for the corresponding NLEE. The effect 
of reductions on the classical r-matrix is discussed. The last Section contains some conclusions. 

2 Preliminaries 

2.1 Cartan— Weyl basis and Weyl group for so(2r) 

Here we fix the notations and the normalization conditions for the Cartan-Weyl generators of 
q ~ so(2r), see e.g. [12]. The root system A of this series of simple Lie algebras consists of the 
roots A = {±(ej — ej),±(ej + ej)} where 1 < i < j < r. We introduce an ordering in A by 
specifying the set of positive roots A + = {ej — ej, ei + ej} for 1 < i < j < r. Obviously all roots 
have the same length equal to 2. 

We introduce the basis in the Cartan subalgebra by h k G f), k = 1, . . . , r where {h k } are the 
Cartan elements dual to the orthonormal basis {e k } in the root space E r . Along with h k we 
introduce also 

r 

H a = y^(a, e k )h k , a G A, 
fc=i 
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where (a, e^) is the scalar product in the root space E r between the root a and e^. The basis 
in so(2r) is completed by adding the Weyl generators E a , a G A. 

The commutation relations for the elements of the Cartan-Weyl basis are given by [12] 



[hk,E a ] — (a,ek)E a , [E a ,E_ a ] — H a , 

' N a ^E a+/3 for a + (3 G A, 
for a + £ A U {0}. 



We will need also the typical 2r-dimensional representation of so(2r). For convenience we choose 
the following definition for the orthogonal algebras and groups 

Xeso(2r) — >X + S X T S = 0, T £ SO(2r) — >S T T S = f, (2.1) 
where by 'hat' we denote the inverse matrix T = T~ x and 

S ^ £(-l) fc+1 (E k - k + E kk ) = ( ° S ° ) , k = 2r + 1 - k. (2.2) 



k=l 



Here and below by Ej k we denote a 2r x 2r matrix with just one non-vanishing and equal to 1 
matrix element at j, k-th position: (Ej k ) mn = 5j m dkn- Obviously Sq = 11. In order to have the 
Cartan generators represented by diagonal matrices we modified the definition of orthogonal 
matrix, see (|2.ip . Using the matrices Ej k defined by equation (|2.2p we get 

hk = E kk - E kk , E ei - ej = Eij - (-1) J+J £^, 
where k = 2r + 1 — k. 



E ei +ei - Ea - (-lf^E-fi, E_ a - E^, 



We will denote by a = Y e k the r-dimensional vector dual to J G f); obviously J = Yl hk- 

k=i fc=i 
If the root a G A + is positive (negative) then (a, a) > ((a, a) < respectively). The normali- 
zation of the basis is determined by 

E- a = Ea, (E- a ,E a )=2, N- a -p = —N a fi. 

The root system A of g is invariant with respect to the Weyl reflections S a ; on the vectors 
y G E r they act as 

_ _ _ 2(a,y) 

JaU = V ~ ~, r", a G A. 

[a, a) 

All Weyl reflections S a form a finite group W s known as the Weyl group. On the root space this 
group is isomorphic to S r ® (Z^) 7 " -1 where S r is the group of permutations of the basic vectors 
ej G W. Each of the Z2 groups acts on E r by changing simultaneously the signs of two of the 
basic vectors ej. 

One may introduce also an action of the Weyl group on the Cartan-Weyl basis, namely [12] 

SaiHp) = A a HpA a l = H Sa /3, 

S a (Ep) = AaEpA" 1 = n a ^Es a i3, n a ^ = ±1. 

The matrices A a are given (up to a factor from the Cartan subgroup) by 

A a = e Ea e- E ~ a e Ea H A , (2.3) 

where Ha is a conveniently chosen element from the Cartan subgroup such that H\ = 1L. The 
formula (|2.3p and the explicit form of the Cartan-Weyl basis in the typical representation will 
be used in calculating the reduction condition following from (|4,16p . 
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2.2 Graded Lie algebras 

One of the important notions in constructing integrable equations and their reductions is the 
one of graded Lie algebra and Kac-Moody algebras [12J. The standard construction is based 
on a finite order automorphism C G Autg, C N = 1. The eigenvalues of C are to k , k = 
0, 1, . . . , N — 1, where u = exp(27ri/A). To each eigenvalue there corresponds a linear subspace 
g( fc ) £ g determined by 

flW = {A: A G 0, C(X) = u fe A}. 

AT-1 , . 

Then = © (J w and the grading condition holds 

k=0 



l (k) ,Q (n) 



C (fc+n) , (2.4) 



where + n is taken modulo A". Thus to each pair {q,C} one can relate an infinite-dimensional 
algebra of Kac-Moody type Qc whose elements are 

A(A) = ]TA fc A fc , A fe G0«. (2.5) 
k 

The series in (I2.5P must contain only finite number of negative (positive) powers of A and 
g(k+N) = g( fc ). This construction is a most natural one for Lax pairs; we see that due to the 
grading condition (12. 4h we can always impose a reduction on L(X) and M(A) such that both 
U(x,t,X) and V(x,t, X) G 0c- In the case of symmetric spaces N = 2 and C is the Cartan 
involution. Then one can choose the Lax operator L in such a way that 

QG0 (1) , JG0 (O) 

as it is the case in Here the subalgebra 0^ consists of all elements of commuting with J. 

r 

The special choice of J = ^ h k taken above allows us to split the set of all positive roots A + 
into two subsets 

A+ = A+ U A+ , A+ = { ei - ej } i<:j , A+ = {a + ej } i<:j . 

Obviously the elements a G have the property a(J) = (a, a) = 2, while the elements j3 G Ajj" 
have the property (3(J) = (P, a) = 0. In this section we outline some of the well known facts 
about the spectral theory of the Lax operators of the type (jl.ip . 



2.3 The scattering problem for L 

Here we briefly outline the basic facts about the direct and the inverse scattering problems 
P31 D31 DSl D51 El DSl D21 EQl ED E21 E3] for the system KHJ for the class of potentials Q(x, t) 
that are smooth enough and fall off to zero fast enough for x — > ±oo for all t. In what follows we 
treat DIH-type symmetric spaces which means that Q(x, t) is an element of the algebra so(2r). 
In the examples below we take r = 4 and ~ so{8). 

The main tool for solving the direct and inverse scattering problems are the Jost solutions 
which are fundamental solutions defined by their asymptotics at x — > ±oo 

lim i/>(x, X)e lXJx = 1, lim <f>(x, X)e iXJx = 1. 

x — ^oo x — y— oo 

Along with the Jost solutions we introduce 



£(x, A) = ip(x, X)e iXJx , <p(x, A) = 4>{x, A)e' 
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which satisfy the following linear integral equations 

£(x, A) = 1 + i f dye- iXJ ^Q(y)ay, A)e iAJ ^), (2.6) 

J oo 

<p(x,\) = 1 + i f dye- iXJ ^Q(yMy,\)e iXJ ^-y\ (2.7) 



These are Volterra type equations which, have solutions providing one can ensure the conver- 
gence of the integrals in the right hand side. For A real the exponential factors in (|2,6p and (|2.7p 
are just oscillating and the convergence is ensured by the fact that Q(x,t) is quickly vanishing 
for x — > oo. 

Remark 1. It is an well known fact that if the potential Q(x, t) G so(2r) then the corresponding 
Jost solutions of equation (jl.ip take values in the corresponding group, i.e. ip(x, A), 4>(x, A) € 
SO(2r). 

The Jost solutions as whole can not be extended for Im A ^ 0. However some of their columns 
can be extended for A £ C+, others - for A 6 C_. More precisely we can write down the Jost 
solutions i/j(x, A) and 4>(x, A) in the following block-matrix form 

V(*,A) = (hT0r,A)),|V + (x,A)» , 0(x,A) = {\4> + (x,\)),\0-(x,\))) , 

where the superscript + and (resp. — ) shows that the corresponding r x r block-matrices allow 
analytic extension for A € C + (resp. A € C_). 

Solving the direct scattering problem means given the potential Q(x) to find the scattering 
matrix T(A). By definition T(A) relates the two Jost solutions 

0(s,A) =V(*,A)T(A), r(A) = (itfy ~a-(\)) (2 ' 8) 

and has compatible block-matrix structure. In what follows we will need also the inverse of the 
scattering matrix 



V(*,A) = </>(x,A)T(A), T(A) = 

where 



c-(A) d-(A) 
-d+(A) c+(A) 



c-(A) = a + (A)(H + /9-/J+)" 1 = (1 + r + r-)- 1 a+(A), 

d-(A) = a+(A)p-(A)(l + pVr 1 = (1 + r + r")- 1 r+(A)a-(A), 

c+(A) = a-(A)(H + = (1 + r-r+)- 1 a-(A), 

d + (A) = a-(A)p+(A)(H + p~p + y l = (1 + r^r + )- 1 r^(A)a+(A). (2.9) 

The diagonal blocks of T(A) and T(A) allow analytic continuation off the real axis, namely 
a + (A), c + (A) are analytic functions of A for A £ C+, while a - (A), c _ (A) are analytic functions 
of A for A € C_. We introduced also /0 =t (A) and t ± (A) the multicomponent generalizations of 
the reflection coefficients (for the scalar case, see [27 ] [17 ] [28]) 

p ± (A) = 6 ± a ± (A) = c ± d ± (A), ^(A) = a ± b =F (A) = rf^A). 

The reflection coefficients do not have analyticity properties and are defined only for A GR. 
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Prom Remark Q] one concludes that T(X) € SO(2r), therefore it must satisfy the second of 
the equations in (|2.1|) . As a result we get the following relations between c 1 * 1 , d^ 1 and a^, b 

c+(A) = s a + ' T (\)s Q , c~(A) = s a~' T (A)s , 

d+(A) = - S ob +,T (A)s , d"(A) = - So b-' T (A)s , (2.10) 
and in addition we have 

p + (A) = -s p + ' T (A)s , P~(A) = -s 0/ o~' T (A)s , 

r+(A) = -s t + > t (\)s , r"(A) = -s r-' T (A) So . (2.11) 

Next we need also the asymptotics of the Jost solutions and the scattering matrix for A — > oo 

lim <j)(x, X)e iXJx = lim ij,(x, X)e iXJx = 1, lim T(A) = 1, 

A^— oo A-*oo A— >oo 

lim a + (A) = lim c~(A) = 1, lim a~(A) = lim c + (A) = 1. 

A— >oo A^oo A— >oo A— >oo 

The inverse to the Jost solutions ip(x, A) and (p(x, A) are solutions to 

i^-$(x,A)(Q(x)-AJ) = 0, (2.12) 
dx 

satisfying the conditions 

lim e" iAJ:z; ^(x, A) = U, lim e~ iXJx (f>(x, A) = 1. (2.13) 



Now it is the collections of rows of i/j(x, A) and (f>(x, A) that possess analytic properties in A 

nx,X )= , 0(x,A)= + , 

V W (a;,A)|/ ^(0 (x,A)| / 

(^(x,A)| = (a^Va^oVf) (*,A), (<^(*,A)| = (s^f,s^t) (x,A). (2.14) 

Just like the Jost solutions, their inverse (|2.14p are solutions to linear equations (|2.12p with 
regular boundary conditions (|2.13|) ; therefore they can have no singularities on the real axis 
AgK. The same holds true also for the scattering matrix T(A) = 4>(x, X)(p(x, A) and its inverse 
T(A) = 4>(x, X)t/)(x, A), i.e. 

a+(A) = $+(x,\)\<l> + (x,\)), a" (A) = (i>~(x, X)\(j)~(x, A)), 

as well as 

c+(A) = (<^ + (x,A)|V + (x,A)), c-(A) = (^(x,A)|V"(x,A)), 

are analytic for A € C± and have no singularities for A 6 R. However they may become 
degenerate (i.e., their determinants may vanish) for some values A^ G C± of A. Below we briefly 
analyze the structure of these degeneracies and show that they are related to discrete spectrum 
of L. 
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3 The fundamental analytic solutions 
and the Riemann— Hilbert problem 

3.1 The fundamental analytic solutions 

The next step is to construct the fundamental analytic solutions (FAS) x A) of (jl.ip . Here 
we slightly modify the definition in [10] to ensure that x ± (x,X) € SO(2r). Thus we define 

X +(x,A) = (\cP + ),^ + c + )) (x,\) = <P(x,\)S + (\) = ^(x,\)T~(\)D + (\) 1 

X ~(x, A) = (l^-fi-), \(j>-)) (x, A) = </>(x, X)S-(X) = ^(x, \)T+(X)D-(\), (3.1) 

where the block-triangular functions S I± (A) and T^(A) are given by 

s-(^{-A- m l)- ^>=(;-T (A, i- p., 

The matrices D ± (X) are block-diagonal and equal 

+ /a+(A) A D - (X) _(c-(X) 

D {X) -{ c+(X))> D {X) -{ o-(A) 

The upper scripts ± here refer to their analyticity properties for A € C±. 

In view of the relations (12.10P it is easy to check that all factors , and take values 
in the group SO{2r). Besides, since 

T(A) = T~{X)D + {X)S + {X) =T+(X)D-(X)S~(X), 

f (A) = S + (X)D + (X)T~ (X) = S-(X)D-(X)f + (X), (3.3) 

we can view the factors S l± , and as generalized Gauss decompositions (see |12j ) of T(A) 
and its inverse. 

The relations between c ± (A), <* (A) and a ± (A), b ± (A) in equation (|2.9p ensure that equa- 
tions (13. 3p become identities. From equations (13. ip . (13. 2p we derive 

X + (*,A) = X -(x,A)G (A), X ~0,A)=x + (*,A)Go(A), (3.4) 

GoW= (t-I1+ T t-t+)' 6«W-( 1+ --"T) < 3 ' 5 > 

valid for A£t. Below we introduce 

X ± (x,X) = X ± (x,X)e iXJx . 

Strictly speaking it is X^(x, A) that allow analytic extension for A 6 C±. They have also another 
nice property, namely their asymptotic behavior for A — > =Loo is given by 

lim X ± (x,X) = II. (3.6) 

A^oo 

Along with X^(x, A) we can use another set of FAS X^(x, X) = X^(x, A)!) 1 * 1 , which also satisfy 
equation (13. 6p due to the fact that 

lim D ± (X) = IL. 

A— >oo 

The analyticity properties of X^(x, A) and X^(x, A) for A € C± along with equation (13.61) are 
crucial for our considerations. 
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Figure 1. The contours j± 



U 7±c 



3.2 The Riemann— Hilbert problem 

The equations (|3.4p and (|3.5p can be written down as 

X + (x,X) = X-(x,X)G(x,X), Aet, 
where 



G{x,X) =e- iXJx G (X)e iXJx . 



Likewise the second pair of FAS satisfy 
X + (x,X) = X-(x,X)G{x,X), AG 

with 



G(x,X) = e~ iXJx G (X)e 



Go (A) 



1 + p-p + p~ 
p+ t 



(3.7) 



(3.8) 



Equation (|3.7|) (resp. equation (|3.8|) ) combined with (|3.6|) is known in the literature [23] as 
a Riemann-Hilbert problem (RHP) with canonical normalization. It is well known that RHP 
with canonical normalization has unique regular solution; the matrix- valued solutions X^(x, A) 
and Xq(x, A) of (|3.7p . (|3.6p is called regular if detX^(x, A) does not vanish for any A G C±. 

Let us now apply the contour-integration method to derive the integral decompositions 
of X^(x, A). To this end we consider the contour integrals 



01 (A) 



and 



1 
2vri 



1 



7+ 



d/i 
p — X 



X+(x,p) 



1 

2vri 



d/i 
p — X 



dp 



X (x,p), 



p — X 2m L_ p — X 

where A € C + and the contours 7± are shown on Fig. [TJ 



X {x,p), 
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Each of these integrals can be evaluated by Cauchy residue theorem. The result for A 6 C + 

are 

3i(A)=X+(*,A)+Y>es X+{x f +VRes X ~^\ (3.9) 

a 2 (A) = X + (x, A) + V Res + V Res . (3.10) 

The discrete sums in the right hand sides of equations (|3.9|) and (|3.10j) naturally provide the 
contribution from the discrete spectrum of L. For the sake of simplicity we assume that L has 
a finite number of simple eigenvalues A^ G C±; for additional details see [10]. Let us clarify 
the above statement. For the 2x2 Zakharov-Shabat problem it is well known that the discrete 
eigenvalues of L are provided by the zeroes of the transmission coefficients a ± (A), which in 
that case are scalar functions. For the more general 2r x 2r Zakharov-Shabat system the 
situation becomes more complex because now a ± (A) are r x r matrices. The discrete eigenval- 
ues A^ now are the points at which a ± (A) become degenerate and their inverse develop pole 
singularities. More precisely, we assume that in the vicinities of A^ a ± (A), c^(X) and their 
inverse a^A), c^(X) have the following decompositions in Laurent series 

o ± (A) = af + (A - Xf)af + ■■■ , c±(A) = cf + (A - Xf)cf + • • • , 
a ± c ± 

where all the leading coefficients , cp, are degenerate matrices such that 

afaf = afaf = 0, cfcf = cfcf = 0. 
In addition we have more relations such as 



**** + "/ 11 • ^+^4 = 11 • 

that are needed to ensure that the identities a ± (A)a ± (A) = 1, c =t (A)c =l= (A) = 1L etc hold true 
for all values of A. 

The assumption that the eigenvalues are simple here means that we have considered only 
first order pole singularities of ^(A) and cf(X). After some additional considerations we find 
that the 'halfs' of the Jost solutions ^^(x, A)} and \(jr t {x, A)) satisfy the following relationships 
for A = Xf 

\^f{x)cf) = ±\4>f(x)T±), \4>f(x)af) = ±\^f(x)pf), 

where 



\i>f{x)) = \^(x,Xf)), \<t>f(x)) = |^(z,A±)>, 

Pj = C 3 ( 'j = h J <l i ' T j = ".< h J = (1 J C j 

and the additional coefficients bf and df are constant r x r nondegenerate matrices which, as 
we shall see below, are also part of the minimal sets of scattering data needed to determine the 
potential Q(x,t). 
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These considerations allow us to calculate explicitly the residues in equations (|3.9p . (|3.10p 
with the result 

X+(x,n) _ (|0),|^ + (x)t+)) X+(x,fi) _ (|V>+(x)p+),|0)) 

Ai=At A* - A AJ - A (U= a+ ^ - A AJ - A 

X-Qr,/i) _ (|^-(x)t-),|0)) _ (|0), \*I>J (x)rr)) 



=A+ A*" A Aj-A M=A + ^ -A AJ -A 



where |0) stands for a collection of r columns whose components are all equal to 0. 

We can also evaluate 0i(A) and 02(A) by integrating along the contours. In integrating along 
the infinite semi-circles of 7± j0 o we use the asymptotic behavior of X (x, A) and X^(x, A) for 
A — ► co. The results are 

4>(x,^ Jx K(x,fi), (3.11) 
ip(x,n)e^ Jx k{x,fi), (3.12) 



11 + 




r°° dfi 


2vri„ 


i-oo At - A 


11 + 




f°° dfi 


2^i. 


l-oo At - A 



K{x,y) = e~ iflJx K (^)e iflJx , K{x,n) = e~ iflJx K (^)e iflJx , 
K(,,\-( K(„\-( 

where in evaluating the integrands we made use of equations (|2.8p . (|2.9p . (|3.7p and (|3.8p . 

Equating the right hand sides of (|3.9p and (|3.11l) . and (13.101) and (I3.12|) we get the following 
integral decomposition for X^(x, A): 

v 



1 f°° du "Xr(x)KiJx) 

X + (x,X) = t + — — ^-X-(x, fJ i)K 1 (x^)+Y / 3 \- t , ^ 
2vri J„ 00 At - A A^ — A 

1 A 00 du *X+(x)K 2j (x) 
2vn /i - A ^ AJ - A 

where Xp(x) = X^(x,\j) and 

AY,!.'-.! = e-^ J - ('' f ) e^'*, A',,+r) = (J' f ) 

Equations ()3. 13j) . (|3.14p can be viewed as a set of singular integral equations which are 
equivalent to the RHP. For the MNLS these were first derived in |25j . 

We end this section by a brief explanation of how the potential Q(x,t) can be recovered 
provided we have solved the RHP and know the solutions X ± (x, A). First we take into account 
that X ± (x,X) satisfy the differential equation 

i— — + Q(x, t)X ± (x, A) - X[J, X ± (x, A)] = (3.15) 
ax 

which must hold true for all A. From equation (|3.6p and also from the integral equations (|3. 13|) . 
(|3.14|) one concludes that X^(x, A) and their inverse Al (x, A) are regular for A — > oo and allow 
asymptotic expansions of the form 

oo oo 

X ± {x, A) = 11 + ^~ s X s (x), Jt+ix, X) = t + Y^ ^ s X s (x). 

8=1 8=1 

Inserting these into equation (|3.15p and taking the limit A — > oo we get 

Q(x,t) = lim A(J- X ± (x,X)JX ± {x,X)}) = [J,Xt(x)]. (3.16) 

A— >oo 
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4 The generalized Fourier transforms 

It is well known that the ISM can be interpreted as a generalized Fourier [ID] transform which 
maps the potential Q(x, t) onto the minimal sets of scattering data Tj. Here we briefly formulate 
these results and in the next Section we will analyze how they are modified under the reduction 
conditions. 

The generalized exponentials are the 'squared solutions' which are determined by the FAS 
and the Cartan-Weyl basis of the corresponding algebra as follows 

*± = X ± {x,\)E a x ± {x,X), =x ± {xA)E- a x ± {x,\), a E A+ 

4.1 Expansion over the 'squared solutions' 

The 'squared solutions' are complete set of functions in the phase space [ID]. This allows one to 
expand any function over the 'squared solutions'. 
Let us introduce the sets of 'squared solutions' 

{*} = {*} c U {*} d , {*} = {*} c U {*} d , 
{*} c = {*+(*, A), ¥Z a (x,\), i<r, Ael}, 

{*} d ^{*+ ;i (:r), *I a;i (x), *I a;i (x)}J =i , 

{*} c = {*l a (x,A), *-(x,A), i<r, A6l}, 

{*} d = *^ a;i (x), *~ ;j (x), ^ ;j (x)} =i , 



where the subscripts 'c' and 'd' refer to the continuous and discrete spectrum of L. The 'squared 
solutions' in bold-face ... are obtained from ... by applying the projector Pqj, i.e. 
*+(x,A) = P J^+(x,A). 

Using the Wronskian relations one can derive the expansions over the 'squared solutions' of 
two important functions. Skipping the calculational details we formulate the results |10| . The 
expansion of Q(x) over the systems {^ ± } and {^^j takes the form 



■ poo 

Q(x) = - dA £ (r+(A)*+(x, A) - r-(A)*I a (x, A)) 

" J— oo 



-oo 
TV 



+ 2 E E ■ r n:j * IK; !.r;) . (4.1) 

i /"°° 

g (x) = __ / d A ^ Gt>i(A)*t a (x,A)-p-(A)*-(z,A)) 

TV 

" 2 E E (Ptu*±ad(?) + P^*^*)) • ( 42 ) 
fc=1 aGA+ 

The next expansion is of ad j 5Q(x) over the systems {^^j and {vf 1 * 1 } 
i f°° 

a dj 1 5Q(x) = — y dA ^ (5r+(A)*+(x,A)+5r-(A)*: Q (x,A)) 

~°° aeA+ 

TV 

+ E E ( 5W «-A x ) - t'w-^)) > ( 4 -3) 

fc=1 aeA+ 
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ad j 1 5Q(x) 



2tt 



dA W(A)*i a (a:,A) + *p-(A)*-(i,A)) 



«eA+ 



where 



N 



(4.4) 



± 



SXfp^j^.jix) + 5p% j ^ a . j (x) 



and = *| a (x, Aj), * ±ay -(x) = d x ®± a (x, A)| A=A ±. 

The expansions (|4.ip . (14. 2p is another way to establish the one-to-one correspondence bet- 
ween Q(x) and each of the minimal sets of scattering data Ti and T2 (|4.6p . Likewise the 
expansions (|4.3p . (14. 4j) establish the one-to-one correspondence between the variation of the 
potential 5Q(x) and the variations of the scattering data 57\ and 572- 

The expansions (|4.3p . (|4.4p have a special particular case when one considers the class of 
variations of Q(x,t) due to the evolution in t. Then 

SQ(x, t) = Q(x, t + St)- Q{x, t) = ^St + (0)((5tf). 

Assuming that 5t is small and keeping only the first order terms in 5t we get the expansions for 
ad j Qt- They are obtained from (14. 3p . (|4.4h by replacing <5p^(A) and <5t^(A) by OtPaW an d 

4.2 The generating operators 

To complete the analogy between the standard Fourier transform and the expansions over the 
'squared solutions' we need the analogs of the operator Dq = — id/dx. The operator Dq is the 
one for which e lXx is an eigenfunction: Dqq iXx = Ae lAx . Therefore it is natural to introduce the 
generating operators A± through 



(A + -A)*t a (x,A) = 0, (Ah 
(A_-A)<&+(x,A) = 0, (A_ 



-A)*"(x,A)=0, 
A)*I Q (x,A) =0, 



(A + 
(A+- 



3 1 ~ =F";j v 



= 0, 

0, 



where the generating operators A± are given by 



A±X(x) = adT 1 i- h i 



dA 



dx 



Q(x), / dy[Q(y),A(2/)] 
±00 



(4.5) 



The rest of the squared solutions are not eigenfunctions of neither A + nor A_ 



(A+ - A+)*I ay .(x) = Vt^ix), (A+ - AJ)* ay .(x 
(A_ - At)*+ W .(x) = *+,.(*), (A_ - AJ)*^(x) = *-,(x), 



i.e., 4f a .j(x) and 3> Q;J (x) are adjoint eigenfunctions of A + and A_. This means that A - , j = 
1,...,N are also the discrete eigenvalues of A± but the corresponding eigenspaces of A± have 
double the dimensions of the ones of L; now they are spanned by both .(x) and 4f^. a .j(x). 
Thus the sets and {<£>} are the complete sets of eigen- and adjoint functions of A + and A_. 
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4.3 The minimal sets of scattering data 

Obviously, given the potential Q(x) one can solve the integral equations for the Jost solutions 
which determine them uniquely. The Jost solutions in turn determine uniquely the scattering 
matrix T(A) and its inverse T(A). But Q(x) contains r(r — 1) independent complex-valued 
functions of x. Thus it is natural to expect that at most r(r — 1) of the coefficients in T(A) 
for A £ R will be independent; the rest must be functions of those. The set of independent 
coefficients of T(A) are known as the minimal set of scattering data. 

The completeness relation for the 'squared solutions' ensure that there is one-to-one corre- 
spondence between the potential Q(x,t) and its expansion coefficients. Thus we may use as 
minimal sets of scattering data the following two sets T« = Tj iC U T^a 

r i N 

7 1 , c = {p + (X),p-(X), Ael}, 7 i,a = {Pj> x f}. =1 > 

7 2 , c = {t + (X),t~(X), A6l}, 7 1)d = [r 3 *A*} , ( 46 ) 

where the reflection coefficients p^(X) and r =t (A) were introduced in equation (|2.9p . A^ are 
(simple) discrete eigenvalues of L and p^ and r- characterize the norming constants of the 
corresponding Jost solutions. 

Remark 2. A consequence of equation PTC]) is the fact that 5 ± (A),T ± (A) G SO(2r). These 
factors can be written also in the form 

S ± (X) = exp I T «W E ±« I » T± W = ex P ( E P^WE± a I • 

\aGA+ / \a£A+ / 

Taking into account that in the typical representation we have E± a E±p = for all roots a, (3 G 
A^ we find that 



E ^( A )^=(o T+ o A) )' £ T « {X)E - 



Ot+(A)^ ^ _ mjJ1 ( s 

,r-(A)0y 

Y. P nx) E±a = (° " + w) , y. = (/ (A) o) • <■") 

a£A+ aGA+ 

where A] 1 " is a subset of the positive roots of so(2r) defined in Subsection 12.21 The formu- 
lae (|4.7p ensure that the number of independent matrix elements of t + (A) and r~(A) (resp., 
p + (X) and p~(X)) equals 2jA^| = r(r — 1) which coincides with the number of independent 
functions of Q(x). 



An important consequence of the expansions is the theorem [10] 

Theorem 1. Any nonlinear evolution equation (NLEE) integrable via the inverse scattering 
method applied to the Lax operator L (jl.ip can be written in the form 



iadj 1 ^ 



where the function /(A) is known as the dispersion law of this NLEE. The generic MMKdV 
equation is a member of this class and is obtained by choosing f(X) = — 4A 3 . // Q(x,t) is 
a solution to (|4.8p then the corresponding scattering matrix satisfy the linear evolution equation 

i^ + f(X)[J,T(X,t)]=0, (4.9) 
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or equivalently 

do ± dXf dp^; , 

[ ~dT T 2fo ^ p = °> -at = °' i_ df T 2/o(A ^ )p S = °' 

dr ± , dA± dr.f , , 

i_ dT ± 2/ °( A ) r = °- "IT = °' i_ df ± 2/o(A ? )r J = °' 

and vice versa. In particular from (|4.9p there follows thata il (X) andc^(X) are time-independent 
and therefore can be considered as generating functionals of integrals of motion for the NLEE. 

Let us, before going into the non-trivial reductions, briefly discuss the Hamiltonian formula- 
tions for the generic (i.e., non-reduced) MMKdV type equations. It is well known (see [10] and 
the numerous references therein) that the class of these equations is generated by the so-called 
recursion operator A = 1/2(A+ + A_) which is defined by equation (|4.5p . 

If no additional reduction is imposed one can write each of the equations in (|4.8p in Hamil- 
tonian form. The corresponding Hamiltonian and symplectic form for the MMKV equation are 
given by 

1 f°° 

^MMKdV = 7 / dx ( tr (JQxQxx) ~ 3tr (JQ' 3 Q X )) , (4.10) 

J — oo 

i f°° / \ i r°° 

n(°) = — / ds tr I ad ^SQix) A [J, ad 7^Q(x)] = — / tr (J<5Q(x) A <5Q(x)). 

1 ./-oo V J ^ J -oo 

The Hamiltonian can be identified as proportional to the fourth coefficient I4 in the asymptotic 
expansion of ^4 + (A) (|5.15p over the negative powers of A 

00 

A + (A) = ^i/ fc A- fc . 

k=l 

This series of integrals of motion is known as the principal one. The first three of these integrals 
take the form 

h = j dxtr(Q 2 (x,t)), I 2 = —\ dxtr(Qad7 1 Q :(; ), 

I 3 = --/ dxtr(QQ xa; + 2Q 4 ), h = — dx (tr (JQ X Q XX ) - 3tr (JQ 3 Q X )) . 
° J -00 31 J— 00 

We will remind also another important result, namely that the gradient of Ik is expressed 
through A as 

V QT{x) I k = --K k - l Q(x,t). 

Then the Hamiltonian equations written through £l( ' and the Hamiltonian vector field X H (o) 
in the form 

n<®(;X Hm )+6H<®=0 (4.11) 

for given by (|4.10p coincides with the MMKdV equation. 

An alternative way to formulate Hamiltonian equations of motion is to introduce along with 
the Hamiltonian the Poisson brackets on the phase space M which is the space of smooth 
functions taking values in g( ' and vanishing fast enough for x — > ±00, see (15. ip . These brackets 
can be introduced by 



{F, G} (0) = i I dxti (V QT(X) F, \j, V qT(x) g\ ) 
J — 00 
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Then the Hamiltonian equations of motions 

^ = fe,tf (0) }( 0) , ^ = te;,# (0) } ( o) (4.12) 

with the above choice for again give the MMKdV equation. 

Along with this standard Hamiltonian formulation there exist a whole hierarchy of them. 
This is a special property of the integrable NLEE. The hierarchy is generated again by the 
recursion operator and has the form 

#S Kd v = -8i4+m, ^ M = - i J°°dx tr (adjHQix) A [J, A m ad ?6Q(x)] ) . 
Of course there is also a hierarchy of Poisson brackets 

J, Ar m V Q T^G 



{F, G} (m) = i I dxtx (v QT(x) F, 
J — 00 ^ 



For a fixed value of m the Poisson bracket {•, -}( m ) is dual to the symplectic form fj( m ) in the 
sense that combined with a given Hamiltonian they produce the same equations of motion. Note 
that since A is an integro-differential operator in general it is not easy to evaluate explicitly its 
negative powers. Using this duality one can avoid the necessity to evaluate negative powers of A. 
Then the analogs of (14. lip and (14.121) take the form: 

rt m \;X Him) ) + 5H^ =0, (4.13) 

^ = { qij ,H(-^} {m) , ^ = K,if(-)} M , (4.14) 

where the hierarchy of Hamiltonians is given by: 

= -4j2fkh+i- m . (4.15) 

k 

The equations (|4.13p and (|4.14p with the Hamiltonian ff(' m ) given by (|4.15p will produce the 
NLEE (|4.8p with dispersion law /(A) = J2 k f k \ k for any value of m. 

Remark 3. It is a separate issue to prove that the hierarchies of symplectic structures and 
Poisson brackets have all the necessary properties. This is done using the spectral decompositions 
of the recursion operators A± which are known also as the expansions over the 'squared solutions' 
of L. We refer the reader to the review papers [261 [TO] where he/she can find the proof of 
the completeness relation for the 'squared solutions' along with the proof that any two of the 
symplectic forms introduced above are compatible. 

4.4 The reduction group of Mikhailov 

The reduction group Gr is a finite group which preserves the Lax representation (|l.ip . i.e. it 
ensures that the reduction constraints are automatically compatible with the evolution. Gr must 
have two realizations: i) Gr C Autg and ii) Gr C Conf C, i.e. as conformal mappings of the 
complex A-plane. To each 5% € Gr we relate a reduction condition for the Lax pair as follows [5] 

C k (L(T k (X))) = Vk L(X), C k (M(T k (X))) = Vk M(X), (4.16) 

where C k € Aut q and Tfc(A) € ConfC are the images of g k and r\ k = 1 or — 1 depending on the 
choice of C k . Since Gr is a finite group then for each g k there exist an integer N k such that 

N k -n 

9k = a - 
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More specifically the automorphisms Ck, k = 1, ...,4 listed above lead to the following 
reductions for the potentials U (x, t, A) and V(x, t, A) of the Lax pair 



U(x,t,X) = Q(x,t) - AJ, V(x,t,X) =^2\ k V k (x,t)-4\ 3 J, 



k=0 



of the Lax representation 

1) C 1 (U\k 1 (X))) = U(X), 



C 1 (V\k 1 (X)))=V(X), 



2) C 2 (U 1 (k 2 {X))) = -U(X), C 2 {V 1 {k 2 {X))) = -V{X), 

3) C 3 (U*( Kl (X))) = -U(X), C 3 (V*(k 1 (\))) = -V(A), 

4) C 4 (U(k 2 (X))) = U(X), C 4 (V( K2 (X)))=V(X), 

where 

a) Kl (X) = X*, b) k 2 (A) = -A. 
The condition (|4.16p is obviously compatible with the group action. 



(4.17) 
(4.18) 
(4.19) 
(4.20) 



5 Finite order reductions of MMKdV equations 

In order that the potential Q(x,t) be relevant for a Dili type symmetric space it must be of 
the form 

Q(x,t) = ^2 (q a (x,t)E a +p a (x,t)E- a ) , 



or, equivalently 

Q(x,t) = ^2 {qijfatfEei+ej ■ +Pij(x,t)E_ ei - ej ) . 

l<i<j<r 

In the next two subsections we display new reductions of the MMKdV equations. 



(5.1) 



5.1 Class A Reductions preserving J 

The class A reductions can be applied also to the MMKdV type equations. The corresponding 
automorphisms C preserve J, i.e. C~ 1 JC = J and are of the form 

C^Wix, \*)C = U(x, A), U(x, A) = Q(x, t) - A J, 

where J is an element of the Cartan subalgebra dual to the vector e± + e 2 + e% + e 4 . In the 
typical representation of so(8) U(x, A) takes the form 



U(x,t,X) 



p(x, t) 



All q(x,t) 
p{x,t) -All J ' 



/ P14 P24, £>34 \ 

P13 P23 p 34 

P12 P23 -P24 

V p± 2 -pi-3 pu J 



q(x,t) 



( qu qi3 qu \ 

924 923 qi 2 

934 q 23 -913 

V g 34 -<?24 914 / 
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Remark 4. The automorphisms that satisfy C 1 JC = J naturally preserve the eigensubspaces 

of adj; in other words their action on the root space maps the subsets of roots A^ onto 

± - A ± 
1 l * 



themselves: CA ± - A± 



We list here several inequivalent reductions of the Zakharov-Shabat system. In the first one 
we choose C = Co to be an element of the Cartan subgroup 



C = exp ^7ri ^2 s khkj 



where take the values and 1. This condition means that Cq = 11, so this will be a Z2- 
reduction, or involution. Then the first example of ^-reduction is 

C - 1 g t (x,t)C = Q(x,t), 

or in components 

Pij — ^ijliji ^ij — ^i^j'j £j — ^ J — il. 

Obviously ej takes values ±1 depending on whether Sj equals or 1. 

The next examples of ^-reduction correspond to several choices of C as elements of the Weyl 
group eventually combined with the Cartan subgroup element Co 

Ci = S ei — e2 S e3 — e4 Co, 

where S ei - ej is the Weyl reflection related to the root ei — ej. Again we have a ^-reduction, 
or an involution 

Cf 1 Q t (x,t)Ci = Q(x,t). 

Written in components it takes the form (ei2 = £34 = 1) 

P12 = ~q*2, P24 = -£23913^ P23 = -£13914, 

P14 = -£13<723> PV3 = ~£23924> P34 = 

The corresponding Hamiltonian and symplectic form take the form 

d x ql 2 d x q 12 + d x q^d x q 34: + £13(^23^914 + d x ql 4 d x q 23 ) 

+ £23 {d x q 24 d x qi 3 + d x ql 3 d x q 24 )j dx + / U12912912 + 934934 

' J— 00 ^ 
\ 2 

+ ^13 (923914 + 924913) + £23(9l4923 + 9l3924)J dx 

/oo 
|<?13<724 + 912934 - 914923 1 2 dx, 
-00 

1 f°° / 

= - / Uq* 12 A 5q 12 + £13(923 A 5q 14 + <5g 14 A <5g 23 ) 
1 J— 00 

+ £23(^24 A S H3 + <5g* 3 A 5q 2A ) + 5q* u A 5q 34 ^jdx. 
Another inequivalent examples of Z2-reduction corresponds to 
C2 = S ei - e2 Co. 
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The involution is C 2 (x,t)C2 = Q(x,t), or in components it takes the form 



P12 = -ei2?12, P24 = -ei3Ql4, P23 = -£l49l3; 

PU = -e23<?24> P13 = -e24<?23: P34 = -£34934- 



As a consequence we get 

8/3 = - / (ei2d x gi 2 d x <?i2 + £34^^34^934 + ei 4 5 x g 13 9 a; g23 
J— 00 

+ £23^924^914 + e24<9x<?235r<?13 + £13^9*4 3e924) <&C 

f ei 2 |gi2 | 2 + £34|g34| 2 + £23924914 + £24923913 + £l49l3923 + £139*4924 J dx 

-00 



-00 

r-oo 



+ £12£34 / 1913924 + 912934 - 9l4923| 2 ^, 
J — 00 

! roo 

= ~ I ( £12^912 A ^912 + £34934 A <$934 + £23^924 A ^914 
1 J —00 

+ £l4^9i 3 A ^923 + £24<5923 A ^913 + £13^9*4 A #924 W- 



Next we consider a Z 3 -reduction generated by C 3 = S ei - e2 S e2 - e3 which also maps J into J. 
It splits each of the sets into two orbits which are 

(O)f = {±(ei + e 2 ) ± (e 2 + e 3 ) ± (ei + e 3 )}, 
(0)f = {±(ei + e 4 ) ± (e 2 + e 4 ) ± (e 3 + e 4 )}. 

In order to be more efficient we make use of the following basis in 
2 2 

where u = exp(2-7ri/3) and a takes values e± + e2 and ei + e 4 . Obviously 

Cg-^WCg = o; fe £W, C^^^Cg = ^i fe) . (5.2) 

In addition, since lo* = u~ l we get (£q )^ = 94°^ and (£a )^ = 3~i 3 ^ for k = 1,2. Then we 
introduce the potential 



QM) = Y,Y1 fe fe) M)4 fe) +A k) (x,t 



k=0 a 

In view of equation (j 5 . 2 j) the reduction condition (|4.17p leads to the following relations between 
the coefficients 

p { 3 = fegr , * = ^r fe) r, = " fc (pi 3 2 - fe) r, 

pS } = ri? = ^g^r , ?a } = ^r fc) r, (5-3) 

where k = 1,2. It is easy to check that from the conditions (|5.3p there follows p\ 2 = 9 12 = 
Pj4 = q$ = 0. So we are left with only one pair of independent functions q^ and qf^ and 
their complex conjugate pf^, q$ . 
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Similarly the reduction (|4.18|) leads to 

p { $ = -(pgr, P ( S = -v- fc (p8r fc) r, *&> = -^(pSt*')*, (5.4) 

where fe = 1,2. Again from the conditions (|5.4p there follows p^ = <7i 2 = P14 = = 0- So 
we are left with two pairs of purely imaginary independent functions: , qyl and pg , ?14 • 

The corresponding Hamiltonian and symplectic form are obtained from the slightly more 
general formulae below by imposing the constraints (|5.3p and (|5.4p . Here for simplicity we skip 
the upper zeroes in qij and jPy 

1 /°° 

-H'MMKdv = g J dx (dlq 12 d x p 12 - d x q 12 (%pi2 + d%qud x pu - d x qud%pu) 
1 f°° 

12 J ( p ii q ii dx ~ Pi^dxQn + QudxPu ~ PndxQn) dx, 

roc g poo 

8/3=0 / dx {d x qi2d x pi2 + d x qud x pu) - - / (g? 4 Pi4 + Q'LpL) dx, 



oc 
00 



= - / dx (Sq 14 A 5pi4 + Sq u A <5pi 2 ) 

i.e. in this case we get two decoupled mKdV equations. 

The Z4-reduction generated by C4 = S ei - e2 S e2 - e3 S e3 - e4 also maps J into J. It splits each 
of the sets into two orbits which are 

(0)f = {±(ei + e 2 ) ± (e 2 + e 3 ) ± (e 3 + e 4 ) ± (ei + e 4 )}, 
(0)± = {±(ei + e 3 ) ± (e 2 + e 4 )}. 

Again we make use of a convenient basis in g(°) 

3 3 

p=0 p=0 

where a takes values e\ + e 2 and ei + e 3 . Obviously 

C 4 - 1 £Wc 4 = i fc £^ ) , C± 1 ?Wc 4 = i k ?W (5.5) 

and in addition, (g^V = ?i 0) , (fi^V = ? ( a~ k) and (£i fe) )* = £^ 4 ~ fc) for k = 1,2,3. Then we 
introduce the potential 

3 



k=0 a 



Q( x ,t) = +p£H*,t 



In view of equation (|5.5p the reduction condition (|4.17p leads to the following relations between 
the coefficients 

p^ = (d 0) r, P {k) = i%t k) r, ^ = i k (pt k) y m 

for k = 1, 2, 3. Here p a , q a coincide with pi 2 , g 12 (resp. pi 3 , qi 3 ) for a = ei+e 2 (resp. a = ei+e 3 ). 
Analogously the reduction (|4.18p gives 

</i 0) = -(^ 0) r, pi 0) = -(pL 0) )*, ?< fc) = -**(?£-*>)•, P {k) = -* k (P { t k) y (5.7) 
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for k = 1,2,3. Both conditions O, lead to p$ = q$ = pff = q[f = for k = 1,3. 

In addition it comes up that £^ = £^ = 9^ = . So we are left with only two pairs of 
independent functions p^, g^ and p^, g^- We provide below slightly more general formulae 
for the corresponding Hamiltonian and symplectic form which are obtained by imposing the 
constraints (|5.6p or (|5.7|) ; again for simplicity of notations we skip the upper zeroes in q^) 

and pfj and replace gj^ and p^ by g*y and pij 

i r 00 

-f^MMKdV = t / dx {dlqi 2 d x pi2 - d x qi 2 dlpi2 + dlq 12 d x p 12 - d x qi 2 d x pu) 

(d x {p\ 2 ) + <9x(j?? 2 )) (g? 2 + g? 2 ) +Pi2Pi25x (gi 2 gi 2 ) 
+ (dx{ql 2 ) + d x (ql 2 )) (p\ 2 +p[ 2 ) + q\ 2 q\ 2 d x {p\ 2 p\ 2 ) J cfc, 

/oo 
{d x qi 2 d x pi 2 + d x q 12 d x p 12 ) dx 
-oo 

- / ((gi2Pi2 + qupu) 2 + (gi2Pi2 + gi2Pi2) 2 ) efo, 

J— oo 

/"OO 

= 2 / (<% 2 A 5pi2 + 5q 12 A Jpi 2 ) . 



3 
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oo 

OO 



Now we get two specially coupled mKdV-type equations. Now we get four specially coupled 
mKdV-type equations given by 

3 3 
dtqo + d x q + -(<?2f>2 + qoPo)d x qo + ~(q 2 po + qop 2 )d x q 2 = 0, 

3 3 
dtq2 + dlq 2 + - (q 2 p + qop 2 )d x q$ + - {q 2 p 2 + qoPo)d x q 2 = 0, 

3 3 
d t po + d x p + -{q 2 p 2 + qoPo)d x p + -{qm + q2Po)d x p 2 = 0, 

3 3 
d t p 2 + d x p 2 + -(q 2 po + qop 2 )d x p + -(qoPo + q2P2)d x p 2 = 0, 

where we use for simplicity 

<7l2 = 90, <Zl2 = 92, Pl2 = P0, Pl2 = P2 

and with second reduction (|5.7p Po = Qo and P2 = — g 2 we have 

3 3 
d t q + d x q + ~(~q2q 2 + qoql)d x qo + -(g 2 go - Qoq 2 )d x q 2 = 0, 

3 3 
d t q 2 + d x q 2 + -(g 2 go - gog 2 )<9xgo + r (<?2<7 2 - qoq 2 )d x q 2 = 0. (5.8) 

Obviously in the system (|5.8p we can put both go, g 2 real with the result 

3 3 
d t q + dlqo +2^1- (&)dxqo = 0, d t q 2 + d z x q 2 + -{q\ - q 2 )d x q 2 = 0. 

The reduction (|4.18p means that gi°^ = igg, p^ 1 = ipo, Qa^ = g 2 , Pa ^ = P 2 with real valued 
Pi ', gV ', i = 0, 2. Thus we get 

dt?o + 3foo + + 1oP2)dxq 2 + ^(p 2 q 2 - goPo ) d xQo = °, 
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« + + \{q^ 2 - q^)d x q\ + + q^)d x q^ = 0, 

n V , <i3 V , ^ / V V V V \ o V i " / VV, V V \ q V n 

OtPo + d x p + -[q 2 Po - q P2)dxPo + ^~1oPo + Q2P2)9xP 2 = 0, 
<9<P2 + <9^2 + ^(Q2P2 - QoPo)dxPo + ^(QoP2 ~ Q2Po) d xP2 = 0- 
5.2 Class B Reductions mapping J into —J 

The class B reductions of the Zakharov-Shabat system change the sign of J, i.e. C" 1 JC = — J; 
therefore we must have also A — ► — A. 

Remark 5. Note that adj has three eigensubspaces W a , a = 0,=L1 corresponding to the 
eigenvalues and ±2. The automorphisms that satisfy C~ l JC = —J naturally preserve the 
eigensubspace Wq but map W_i onto Wi and vice versa. In other words their action on the 
root space maps the subset of roots Af onto and vice versa: CAf = Af. 

Here we first consider 

Obviously the product of the the above 4 Weyl reflections will change the sign of J. Its effect 
on Q in components reads 

* * _ 

Qij ^ijQijj Pij ^ijPiji ^ij ^i^j: 

i.e. some of the components of Q become purely imaginary, others may become real depending 
on the choice of the signs €j. 

There is no Z3 reduction for the so(8) MMKdV that maps J into —J. So we go directly 
to the Z4-reduction generated by Cq = S ei+e2 S e2+e3 S e3+e4 which maps J into — J. The orbits 
of C3 are 

(0)f = {±(ei + e 2 ) =F (e 2 + e 3 ) ± (e 3 + e 4 ) =F (ei + e 4 )}, 
(0)± = {±(ei + e 3 ) T (e 2 + e 4 )}- 

Again we make use of a convenient basis in 

3 3 

p=0 p=0 

where a takes values e\ + e 2 and e\ + e 3 . Obviously 

C^llWCe = i fc £( fc) , C^Ce = i fc J«, (5.9) 

where again = ^ , (fi^V = 3 { a~ k) and (8^)* = £^" fc) for A; = 1,2,3. Then we 

introduce the potential 



Q(x,t) =^^(gW(,,t)£W +p W(M 



fc=0 a 



In view of equation (|5.9p the reduction condition (|4.17p leads to the following relations between 
the coefficients 



P 



L 0) = (<z£ 0) r, *4 fc) = g£ fc) =**(?#-*>)*, ^ = 1,2,3, (5.10) 
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while the reduction (|4.18p gives 



(0) = -(9i 0) r, « { a k) = -i k ^- k) r, fc = i,2,3, 



(0) --(pL 0) r, P {k) = ~i k ¥t k) )\ k = 1,2,3, (5.11) 



where a takes values e\ + e2 and e\ + e3. From the conditions (|5.10p there follows p k ^ = = 
Pi3^ = 9i3 = for /c = 1,3. In addition however, it comes up that £^g = £^ = ujg = 
So we are left with only two pairs of independent functions pf 2 , q^ and p^ , gj^- We provide 
below slightly more general formulae for the corresponding Hamiltonian and symplectic form 
which are obtained by imposing the constraints (|5.1U|) or (|5.11|) and again for simplicity we skip 

the upper zeroes in q^ and pij and replace q^ and pf 2 by q\ 2 and p\ 2 
1 f°° 

-f^MMKdv = t / ete (dlq 12 d x pi2 - d x q\ 2 d 2 x p\ 2 + d 2 x q\ 2 d x px 2 - d x qi 2 d x p 12 ) 
./ — oo 

(^(pL) + d x(qj 2 )) {qj 2 +Pi 2 ) + d x (qnpn) P\ 2 q\ 2 
~ {px{q\ 2 ) + d x (p\ 2 )) (p\ 2 + gf 2 ) - d x (pi 2 qi 2 ) qi 2 pi 2 ^jdx, 

/oo 
(d x qi 2 d x pi 2 + d x q X2 d x p\ 2 ) dx 
-oo 

- / ( (<?i2Pi2 + qupuf + (qi2Pi2 + qupu) 2 ) dx, 

/oo 
dx (5q 12 A <5pi2 + Sq 12 A <5pi 
-oo 



_3_ 
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oo 

oo 



>12) ■ 

Now we get two specially coupled mKdV-type equations. 

The class B reductions also render all the symplectic forms and Hamiltonians in the hierarchy 
real-valued. They allow to render the corresponding systems of MMKdV equations into ones 
involving only real-valued fields. 'Half of the Hamiltonian structures do not survive these 
reductions and become degenerate. This holds true for all symplectic forms fi( 2m ) and integrals 
of motion I 2m with even indices. However the other 'half of the hierarchy with Q( 2m+1 ) and 
integrals of motion I 2m +i remains and provides Hamiltonian properties of the MMKdV. 

Now we get four specially coupled mKdV-type equations given by 

3 3 
d t qo + d x q + -(<7o<72 + PoP2)d x p 2 + ~{p 2 q 2 + qoPo)d x qo = 0, 

3 3 
d t q 2 + dlq 2 + - (q 2 p 2 + qoPo)d x q 2 + -(<?o<?2 + PoP2)d x p = 0, 

3 3 
d t po + dlpo + -{q 2 p 2 + qoPo)d x p + -(<?o<?2 + PoP2)d x q 2 = 0, 

3 3 
d t p 2 + dlp 2 + -{q 2 q +PoP2)d x q + -(qoPo + q2P2)d x p 2 = 0, 

where we use for simplicity 

qi2 = qo, qi2 = <?2, pu = po, P12 = P2 

and with second reduction (|5.10p po = q$ and p 2 = —q 2 we have 

3 3 
dtqo + d x q - -(<?o<?2 - qlq 2 )d x ql - -{q 2 q 2 - qml)d x qo = 0, 

3 3 
d t q 2 + dlq 2 - -(q 2 q 2 - qoql)d x q 2 + -(q q2 - qlql)d x q* Q = 0. (5.12) 
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Obviously in the system (|5.12p we can put both qo, q 2 real with the result 

d t q + dlq + -{ql - q 2 )d x q = 0, d t q 2 + d x q 2 + ^(tfo ~ 92)^r?2 = 0. 

The reduction (|5.1ip means that q^ = iq^ , Pa = iPo , Qa^ = q 2 -> Pa = P 2 with real valued 
p( , q( , i = 0, 2. Thus we get 

« + d 3 x q V o + l(QoQ V 2 + Vlpi)d x pi + \(pW 2 ~ «)^ V = 0, 
dtqi + fl££ + ^(g 2 >v _ g VpV Rg y _ 3^ +p Y p v ) ^v = Qj 
+ ^ V + |(«#pV _ ^v^v + 3 (? v g v + p v p v } ^ g v = Q) 

a V , a3 V 3 - v V i V V \ o V 3, V V V V\n V n 
<9tP 2 + ^2 - 2^2 90 +P0P2)dxq ~ ^{QoPO -Q2P2)OxP2 = 0- 

The class B reductions also render all the symplectic forms and Hamiltonians in the hierarchy 
real-valued. They allow to render the corresponding systems of MMKdV equations into ones 
involving only real-valued fields. 'Half of the Hamiltonian structures do not survive these 
reductions and become degenerate. This holds true for all symplectic forms Q^ 2m ^ and integrals 
of motion I 2m with even indices. However the other 'half of the hierarchy with f2( 2m+1 ) and 
integrals of motion Xgm+i remains and provides Hamiltonian properties of the MMKdV. 

5.3 Effects of class A reductions on the scattering data 

The reflection coefficients p ± (A) and T :t (A) are defined only on the real A-axis, while the diagonal 
blocks ^(A) and c ± (A) (or, equivalently, D^(X)) allow analytic extensions for A £ C±. From 
the equations (|2.9p there follows that 

a + (A)cT(A) = (U + / 9-p + (A)r\ a-(A)c+(A) = (1 + pVW)" 1 , (5.13) 
c-(A)a+(A) = (l + r + r-(A))~ 1 , c+(A)a~(A) = (U + r-r+(A))- 1 . (5.14) 

Given 7± (resp. T2) we determine the right hand sides of (I5.13P (resp. (|5.14p ) for A 6 R. 
Combined with the facts about the limits 

lim a + (A) = 11, lim c~(A) = U, lim o~(A) = 1, lim c + (A) = 1, (5.15) 

A— »oo A— >oo A— >oo A— *oo 

each of the relations (|5.13|) . (|5.14|) can be viewed as a RHP with canonical normalization. Such 
RHP can be solved explicitly in the one-component case (provided we know the locations of 
their zeroes) by using the Plemelj-Sokhotsky formulae [23]. These zeroes are in fact the discrete 
eigenvalues of L. One possibility to make use of these facts is to take log of the determinants of 
both sides of (|5.13p getting 

A + (A) + C~(A) = -lndet(l + p~p + (X), A G R, 

where 

A ± (A) =lndeta ± (A), C±(A) =lndetc ± (A). 
Then Plemelj-Sokhotsky formulae allows us to recover A^A) and C ± (A) 

/oo j N \ ^+ 

— ^ In det(l + p-p + (v)) + £ In —-L, (5.16) 
-00 V ~ A ^ A - Xj 
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where A(X) = A + (X) for A G C+ and A(X) = -C~(X) for A € C_. In deriving (pUB]) we have 
also assumed that A^ are simple zeroes of A ± (X) and C ± (A). 

Let us consider a reduction condition (|4.17p with C\ from the Cartan subgroup: C\ = 
diag (-B+, where the diagonal matrices -B± are such that .B± = 1. Then we get the following 
constraints on the sets Ti 2 

p~(X) = (5_p + (A)i? + )t, pj = (S_p+fl+)t, AT = (At)*, 
r-(A) = ( J B + r+(A) J B_)t, r r = (B + r/S_)t, Aj = (A+)*, 

where j = 1, . . . , N. For more details see Subsections 14.31 and 14.41 

Remark 6. For certain reductions such as, e.g. Q = —Q^ the generalized Zakharov-Shabat 
system L{X)ip = can be written down as an eigenvalue problem Htfj = Xtp(x, A) where £ 
is a self-adjoint operator. The continuous spectrum of XL fills up the whole real A-axis thus 
'leaving no space' for discrete eigenvalues. Such Lax operators have no discrete spectrum and 
the corresponding MNLS or MMKdV equations do not have soliton solutions. 

From the general theory of RHP [21] one may conclude that f)5 . 13|) , (|5.14p allow unique solu- 
tions provided the number and types of the zeroes A^ are properly chosen. Thus we can outline 
a procedure which allows one to reconstruct not only T(A) and T(X) and the corresponding 
potential Q(x) from each of the sets Tj, i = 1,2: 

i) Given T2 (resp. Ti) solve the RHP (15. 131) (resp. (|5.14p ) and construct a ± (A) and c ± (A) for 

A G C±. 

ii) Given Ti we determine & ± (A) and ci ± (A) as 

6±(A) = p ± (A)a ± (A), d±(A) = c ± (A)p ± (A), 
or if T2 is known then 

b ± (X) = a ± (A)r ± (A), d ± (A) = r ± (A)c ± (A). 

iii) The potential Q(x) can be recovered from 7\ by solving the RHP (|3.7p and using equa- 

tion (I3T6D . 

Another method for reconstructing Q(x) from Tj uses the interpretation of the ISM as gene- 
ralized Fourier transform, see [27 \ \28\ [29] . 

Let in this subsection all automorphisms C% are of class A. Therefore acting on the root space 

r 

they preserve the vector ^ which is dual to J, and as a consequence, the corresponding Weyl 

k=l 

group elements map the subset of roots A± onto itself. 

Remark 7. An important consequence of this is that C, will map block-upper-triangular (resp. 
block-lower-triangular) matrices like in equation (|3.2p into matrices with the same block struc- 
ture. The block-diagonal matrices will be mapped again into block-diagonal ones. 

From the reduction conditions (|4.17|) - (|4.20p one gets, in the limit x — > 00 that 
a) Ci((«i(A)J)t) = AJ, b) C 2 ((k 2 (A)J) t ) = -AJ, 

c) C 3 ((k 3 (X)JT) = AJ, d) C 4 ((« 4 (A) J)) = AJ. (5.17) 
Using equation (I5.17P and Q(J) = J one finds that: 

a) «i(A) = A*, b) K2 (A) = -A, c) k 3 (A) = -A*, d) /54(A) = A. 
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and 



It remains to take into account that the reductions (|4,17p ~ (|4.20p for the potentials of L lead 
to the following constraints on the scattering matrix T(A) 

a) d(rt(A*)) = f(A), b) C 2 (T T (-X))=f(X), 

c) C 3 ((T*(-A*))=T(A), d) C 4 (T(A))=T(A). 

These results along with Remark [7] lead to the following results for the generalized Gauss 
factors of T(A) 

a) C 1 (S + >\\*)) = S-(\), C 1 (T-'t(A*))=f+(A), 

b) C 2 (S + > T (-\)) = S-(\), C 2 (T- T (-A)) = T+(A), 

c) CsiS^i-X*)) = S ± (X), C 3 (T±>*(-\*)) = T±(A), 

d) C A (S ± (X)) = S ± (\), C 4 (T ± (A))=T ± (A) 

a) C 1 (D + >\\*)) = D-(\), b) C 2 (D + > T (-\)) = D~(\), 
c) C 3 (D ± >*(-\)) = D ± (\), d) C A (D ± (X)) = D ± (\). 

5.4 Effects of class B reductions on the scattering data 

In this subsection all automorphisms Cj are of class B. Therefore acting on the root space they 

r r 

map the vector ^ dual to J into — Yl e k- As a consequence, the corresponding Weyl group 

fc=i k=l 
elements map the subset of roots Af onto Aj~ = — Af. 

Remark 8. An important consequence of this is that C{ will map block-upper-triangular into 
block-lower-triangular matrices like in equation (|3.2p and vice versa. The block-diagonal matrices 
will be mapped again into block-diagonal ones. 

Now equation (|5.17p with Cj(J) = — J leads to 

a) Kl (A) = -A*, b) K2 (A) = A, c) k 3 (A) = A*, d) k 4 (A) = -A. 

The reductions (|4.17p - (|4.20p for the potentials of L lead to the following constraints on the 
scattering matrix T(A): 

a) d(Tt(-A*))=f(A), b) C 2 (T T (X)) = f(X), 
c) C 3 ((T*(A*))=T(A), d) C 4 (r(-A))=T(A). 

Then along with Remark [8] we find the following results for the generalized Gauss factors 
of T(A) 

a) C 1 (5 ± 't(-A*)) = 5 ± (A), C^T^-A*)) = T±(\), 

b) C 2 (5 ± ' T (A)) = 5 ± (A), C 2 (T ± ' T (A))=T ± (A), 

c) C 3 (S + '*(\*)) = S-(\), C 3 (T->*(A*))=T+(A), 

d) C 4 (-S+(-A)) = 5-(A), C 4 (T"(-A)) =T+(A) 



and 



a) C 1 (D ± '\-\*)) = D±(\), b) C 2 (D ± - T (A)) = J D±(A), 
c) C 3 ( J D+'*(A*)) = J D-(A), d) C 4 (Z)+(-A)) = Z)-(A). 
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6 The classical r-matrix and the NLEE of MMKdV type 

One of the definitions of the classical r-matrix is based on the Lax representation for the cor- 
responding NLEE. We will start from this definition, but before to state it will introduce the 
following notation 



[u(x,\)®U(y,n)}, 



which is an abbreviated record for the Poisson bracket between all matrix elements of U(x, A) 
and U(y,fi) 



\u(x,\)®U(y,Li)\ = {U ik (x,X),U lm (y,Li)}. 

y- ' > ik.lm 



ik,lm 

In particular, if U(x, A) is of the form 



U(x, A) = Q(x,t) - A J, Q(x,t) = ^(qirEir +p ri E ri ), 



i<r 



and the matrix elements of Q(x,t) satisfy canonical Poisson brackets 
UksiPri} = i5 ik 5 rs 5(x - y), 

then 

® E ri - E ri (g) E ir )5(x - y). 

i<r 

The classical r-matrix can be defined through the relation [19] 

{U(x, A) ® U(y, M)} = i [r(A - fi), U(x, A) ® U + 1 ® U{y, //)] S(x - y), (6.1) 

which can be understood as a system of iV 2 equation for the N 2 matrix elements of r(A — fi). 
However, these relations must hold identically with respect to A and fx, i.e., (|6.ip is an overdeter- 
mined system of algebraic equations for the matrix elements of r. It is far from obvious whether 
such r(A — /i) exists, still less obvious is that it depends only on the difference A — fi. In other 
words far from any choice for U(x, A) and for the Poisson brackets between its matrix elements 
allow r-matrix description. Our system (|6.ip allows an r-matrix given by 

where P is a constant iV 2 x iV 2 matrix 

N r 

P = ^2 (E a ® E_ a + £_ Q ® E a ) + ^2 h k ® h k . 

«GA+ k=l 

The matrix P possesses the following special properties 
[P,X® 1 + -0.&X] =0 VIeg. 
By using these properties of P we are getting 

[P, Q(x) ® 11 + 11 ® Q{x)\ = 0, (6.3) 
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i.e., the r.h.s. of (|6.ip does not contain Q(x,t). Besides: 
[P, X J ® & + /ill ® J] = (A - n) [P, J ® a] 

= -2(A - /i) ( (E a ® £_ a - £_ a ® £ a ) , (6.4) 

where we used the commutation relations between the elements of the Cartan-Weyl basis. The 
comparison between (|6.3p . (|6.4p and (|6.ip leads us to the result, that r(A — fj,) (|6.2p indeed 
satisfies the definition (16.11). 



Remark 9. It is easy to prove that equation (|6.ip is invariant under the reduction condi- 
tions (|4.17p and (|4.18p of class A, so the corresponding reduced equations have the same r- 
matrix. However (|6.ip is not invariant under the class B reductions and the corresponding 
reduced NLEE do not allow classical r-matrix defined through equation 



Let us now show, that the classical r-matrix is a very effective tool for calculating the Poisson 
brackets between the matrix elements of T(A). It will be more convenient here to consider 
periodic boundary conditions on the interval [— L, L], i.e. Q(x — L) = Q(x + L) and to use the 
fundamental solution T(x, y, A) defined by 

. dT(x y, X) + ^ = Q ^ T (x,x,A) = U. 

dx 

Skipping the details we just formulate the following relation for the Poisson brackets between 
the matrix elements of T(x, y, A) [19] 



|r(x, y, A) <g> T(x, y, //)} = [r(A - fJ,),T(x, y, A) T(x, y, fjt)] . (6.5) 

The corresponding monodromy matrix Tl(X) describes the transition from —L to L and 
Tl(X) = T(—L, L, A). The Poisson brackets between the matrix elements of ?l(A) follow directly 
from equation (|6.5p and are given by 



{T L (X)®T L (v)} = [r(A - fi), T L (X) (8> T L (n)] . 

An elementary consequence of this result is the involutivity of the integrals of motion 1^^ 
from the principal series which are from the expansions of 

oo oo 

lndeta+(A) = ^/ L ,fcA~ fc , -lndetcT^A) = ^/ L>fc A~ fc , 

k=l k=l 
oo oo 

lndetc+(A) = JL,kX~ h , - lndet a~(A) = J L ,kX~ k . (6.6) 
fe=i fe=i 

An important property of the integrals II & and ^ is their locality, i.e. their densities depend 
only on Q and its x-derivatives. 

The simplest consequence of the relation (|6.5|) is the involutivity of Ii.k, Jl,U- Indeed, taking 
the trace of both sides of (|6.5p shows that {tr Tl(X), tr !/,(//)} = 0. We can also multiply both 
sides of (|6.5p by C <g> C and then take the trace using equation (|6.3|) ; this proves 



{trT L (A)C,trT L ( At )C} = 0. 
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In particular, for C = IL + J and C = 11 — J we get the involutivity of 

{tra+(A),tra7»} = 0, {tra7(A),tra7»} = 0, 
{trc+(A),trc+(^)} = 0, {tr c£(X), tr c7»} = 0. 

Equation (|6,5|) was derived for the typical representation of (25 ~ SO(2r), but it holds 
true also for any other finite-dimensional representation of <5. Let us denote by ~ A^V^ 1 -* 
the fc-th fundamental representation of (5; then the element ?l(A) will be represented in 
by A fe T L (A) - the fc-th wedge power of Tl(X), see [12] . In particular, if we consider equation (|6.5p 
in the representation and sandwich it between the highest and lowest weight vectors in 
we get [30] 



{deta+(A),deta+(^)} = 0, {det c7(A), det c7(^)} = 0. (6.7) 

Likewise considering (|6,5p in the representation y( m ) and sandwich it between the highest 
and lowest weight vectors in y( m ) we get 

{deta7(A),deta~(^)} = 0, {det c+(A), det c[{ji)} = 0. (6.8) 

Since equations (|6.7p and (|6.8p hold true for all values of A and we can insert into them the 
expansions (|6.6|) with the result 



{h,kJL, P } = 0, { J L ,fc, Jl, p } = 0, fc,p= 1,2, .... 

Somewhat more general analysis along this lines allows one to see that only the eigenvalues 
of a |~ (A) and c^(A) produce integrals of motion in involution. 

Taking the limit L — > co we are able to transfer these results also for the case of potentials 
with zero boundary conditions. Indeed, let us multiply (|6.5p by E(y, A) <g> E(y, //) on the right 
and by E~ 1 {x, A) <S> E~ l (x,[i) on the left, where E(x,X) = exp(— iAJx) and take the limit for 
x — > oo, y — > — oo. Since 

g ia;(A— /i) 

lim — = ±mS(X — jit), 

x^±oo X — fl 

we get 

|r(A) ® T(/i)} = r+(A - /u)T(A) <g> T(/i) - T(A) ® T(^)r_(A - /i), 

r±(A - /u) = - 2 / A "'_ x f it, hk + 2 ^ + K - Q ® ^) ] ± ivr5(A - /u)n j, 



, A-=l 



where IIoj is defined by 

n 0J = ^(E a ® E_ a - E_ a ® E a ). 

QGA+ 

Analogously we prove that i) the integrals 1^ = lim 1^ & and J p = lim Jl jP are in involution, 



i.e. 



{-ffc, Ip} — {h, Jp} — {Jk, Jp} — 0, 

for all positive values of k and p; ii) only the eigenvalues of a^A) and c^(X) produce integrals 
of motion in involution. 
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7 Conclusions 

We showed that the interpretation of the ISM as a generalized Fourier transform holds true 
also for the generalized Zakharov-Shabat systems related to the symmetric spaces Dili. The 
expansions over the 'squared solutions' are natural tool to derive the fundamental properties 
not only of the MNLS type equations, but also of the NLEE with generic dispersion laws, in 
particular MMKdV equations. Some of these equations, besides the intriguing properties as 
dynamical systems allowing for boomerons, trappons etc., may also have interesting physical 
applications. 

Another interesting area for further investigations is to study and classify the reductions of 
these NLEE. For results along this line for the MNLS equations see the reports |10j and [9]; 
reductions of other types of NLEE have been considered in [261 El El [6]. 

One can also treat generalized Zakharov-Shabat systems related to other symmetric spaces. 
The expansions over the 'squared solutions' can be closely related to the graded Lie algebras, 
and to the reduction group and provide an effective tool to derive and analyze new soliton 
equations. For more details and further reading see [111 126]. 

In conclusion, we have considered the reduced multicomponent MKdV equations associated 
with the Dill-type symmetric spaces using A.V. Mikhailov reduction group. Several examples 
of such nontrivial reductions leading to new MMKdV systems related to the so(8) Lie algebra 
are given. In particular we provide examples with reduction groups isomorphic to Z2, Z3, Z4 
and derive their effects on the scattering matrix, the minimal sets of scattering data and on the 
hierarchy of Hamiltonian structures. These results can be generalized also for other types of 
symmetric spaces. 
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